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ABSTRACT 

We present observations of a new low-mass double-lined eclipsing binary sys- 
tem discovered using repeat observations of the celestial equator from the Sloan 
Digital Sky Survey II. Using near-infrared photometry and optical spectroscopy 
we have measured the properties of this short-period [P=0. 407037(14) d] system 
and its two components. We find the following parameters for the two compo- 
nents: M x = 0.272 ± O.O2OM , R4 = 0.268 ± O.OlOi? , M 2 = 0.240 ± O.O22M , 
R 2 = 0.248 ± 0.0090^, T x = 3320 ± 130K, T 2 = 3300 ± 130K. The masses and 
radii of the two components of this system agree well with theoretical expec- 
tations based on models of low-mass stars, within the admittedly large errors. 
Future synoptic surveys like Pan-STARRS and LSST will produce a wealth of 
information about low-mass eclipsing systems and should make it possible, with 
an increased reliance on follow-up observations, to detect many systems with low- 
mass and sub-stellar companions. With the large numbers of objects for which 
these surveys will produce high-quality photometry, we suggest that it becomes 
possible to identify such systems even with sparse time sampling and a relatively 
small number of individual observations. 
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Introduction 



Studies of detached, double-lined eclipsing binary systems (DEBs) allow for direct mea- 
surements of the masses and radii of stars. These measurements have direct implications for 
theories of stellar structure and evolution. Our understanding of the structure of stars on 
the lower main sequence has improved over the past decade, but relatively few direct mea- 
surements of stellar mass and radius exist for objects with masses O.1M < M* < O.4M . 
Measurements of stars in this mass range allow us to test models of stellar structure where 
they are least well constrained by observation. There are four known main sequence DEBs 
inhabiting this mass range and the properties of these system s are well-measured ( jRibas 
20031 . iMetcalfe et al.lll996l iHebb et al.ll2006l IVaccaro et al.ll2007l ). In order to refine physical 
models of low-mass stars, more systems must be discovered and studied. Large-area synop- 
tic surveys offer the best hope for the identification of new low-mass eclipsing systems. The 
photometric databases of shallow, wide-angle, surveys such as those aimed at the discovery 
of GRB afterglo ws or transiting planets, are already being combed for the most interesting 
binary systems (ILopez-Morales et al.ll2006l ; iDevor &: Charbonneaull2006l ). but the number of 
M dwarfs in these shallow surveys is limited. These surveys have also been used to identify 
low-mass si ngle-lined systems f rom which the mass and radius of an unseen companion can 
be inferred (IBeatty et al.l 120071 ). 



The advent of deeper surv eys with sparser tim e sampling (less than o ne observation 
per day), such as Pan-STARRS (IHodapp et al.ll2004l ) and LSST (jTysonll2002l ). will become a 
vital resource for studies of variable objects of all types. Since low-mass stars are intrinsically 
faint, deep surveys may prove particularly powerful for studies of these objects. The Sloan 
Digital Sky Survey (SDSS) and its follow-up survey SDSS-II are excellent prototypes of 
the types of data that these future surveys will produce. Here, we report the first results 
from a program to identify DEBs in SDSS-II observations of an area of the sky scanned 
repeatedly as part of the SDSS Supernova Survey program. We present the discovery of 
a new DEB along with mass and radius measurements of these two low-mass stars based 
on spectra taken with Keck/LRIS and near infrared (NIR) photometry gathered with the 
robotic telescope PAIRITEL. This discovery helps to fill in a poorly sampled portion of the 
Mass-Radius relation for low-mass stars. In Section 2, we describe the identification of this 
DEB and the follow-up observations. In Section 3, we describe our analysis, modeling, and 
the derivation of the physical parameters of the system. In Sections 4 and 5 we discuss our 
results in the context of the known low-mass eclipsing binary systems and future prospects 
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for the discovery of similar systems. 



Observations 



2.1. Identification of SDSS-MEB-1 



The SDSS-II Supernova Survey (jFrieman et al.ll2008l ) repeatedly observes about 300 deg 2 
of sky centered around the celestial equator. The raw data, both imaging and uncalibrated 
outputs from the photometric pipeline, are made publicly available through the SDSS SN 
Archive ServeiQ. We downloaded all of the imaging data from i band and the corresponding 
pipeline outputs available as of April, 2006. In total, nearly 2 x 10 5 individual 2048 by 1489 
pixel images from 70 observing runs were available. We calibrate d the uncalibrated outputs 



of th e photometric pipeline against stars from the SDSS DR4 (lAdelman-McCarthy et al. 



20061 ) photometric catalog. This was done by matching the pixel coordinates of stars in 
the raw i-band images to the coordinates of preselected catalog F-G-K stars and calculating 
a simple photometric solution for each individual image using raw PSF magnitudes (no 
color terms). Once the positions of stars in the raw i-band images were matched to the 
photometric pipeline out puts, a simil a r pho tome tric solution was calculated in the r— and 
z-bands. Recent work by llvezic et al.l (120071 ) and iPadmanabhan et al.l (120081 ) presents more 
sophisticated methods for calibrating repeat observations in SDSS. We found our technique to 
produce individual light curves with RMS of between 0.03 and 0.05 mag (i <19.0), adequate 
precision for the detection of large eclipse signals. 

We preselected a set of 19,000 M dwarf targets from t he photomet r ic dat abase of SDSS 
DR4. This was done following M dwarf colors presented by lWest et al.l (120051 ) and requiring 
that objects have flags corresponding to high photometric quality. Our targets all have 
i — z > 0.84 and r < 21.2, corresponding roughly to a sample of objects of spectral types 
M4 and later with r band photometric precision better than 5%. Faint objects with colors 
of i — z > 0.8 are expected to be almost exclusively M dwarfs since other types of point 
sources, s uch as qua sars and M giants, are very rare in this portion of the color-magnitude 
diagram (jFanl Il999l ). We generated r,i,z light curves for all of these objects, each light 
curve consisting of between 10 and 30 observations in each band. We identified potential 
eclipsing sources by searching for objects that exhibited decreases in brightness of more than 
0.20 mag which were found to occur simultaneously in all three bands. This conservative 
requirement for the large amplitude of the dimming and the requirement that it be found 
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in all three bands resulted in a robust method for identifying eclipsing candidates. Since we 
have relatively few observations of each target, we require a conservative detection threshold 
to ensure a very low rate of false positives. Our analysis produced a list of 13 candidates 
that have been observed as part of an ongoing program to identify DEBS in SDSS-II. Of 
these candidates, SDSS031824-010018 was the most promising. The SDSS-II light curve of 
this object, which we call SDSS-MEB-1, contained 15 points in r,i, and z and displayed two 
epochs with fluxes that were lower than the mean flux by Am > 0.3 mag in each band. An 
expanded SDSS-II light curve, containing additional data, is shown in Figure [TJ The basic 
observed parameters of SDSS-MEB-1 are given in Table 1. 



2.2. Follow-up Photometry 



We observed SDSS-MEB-1 with the PAIRITEL robotic NIR observatory between 9 
September, 2006 and 8 October, 2006. Observations consisted of 30 to 60 minutes of in- 
dividual 7.8s integrations on each clear night. The PAIRITEL camera takes images in the 
J, H, and K s bands simultaneously. The PAIRITEL imag er was originally used for the 
southern portion of the 2MASS survey (ISkrutskie et al.l 120061 ) . The observations were auto- 



m atically scheduled, c ollec ted, and reduced w ith the fully robotic observing system described 



in 



Bloom et al.l (120051 ) and lBlake et al.l (120051 ). The images were gathered using a randomized 



dither pattern and individual integrations were combined into mosaics with total integration 
time of ~300s using SWarp (IBertin et al.ll2002l ). We produced differential aperture photom- 
etry relative to a set of comparison stars from the 2MASS catalog and photometric errors 
were estimated based on read noise and photon noise. Our analysis resulted in 937 individual 
measurements in J, H, and K s bands. We report the J, H, and K s measurements in Table 
2. 



2.3. Spectroscopy 



We observed SDSS-MEB-1 spectroscopically using the LRIS (jOke et al.lll995l ) instru- 
ment on Keck on 20 September 2006 and 21 November 2006. Spectra were taken at reso- 
lutions R~ 900 and R~ 3500 using the 400 grooves mm -1 , 8500 A blaze and 1200 grooves 
mm -1 , 7500 A blaze gratings with a 1" slit. These two setups cover the spectral ranges 5550- 
9270 A and 6335-7600 A , respectively. We observed the standard star HD 19445 for flux 
calibration. We ex t racte d the red channel spectra following the optimal extraction method 
outlined by iHornd (119861 ). We did not observe any radial velocity standards, so wavelength 
solutions were determined using both HgNeAr comparison lamps and isolated night sky 
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emission lines from the catalog of lOsterbrock et al.l (119971 ). All spectra showed prominent 
emission due to Ha and the R ~ 3500 spectra revealed two distinct emission lines. The 
detection of Ha emission demonstrates that the stars have chromospheric activity. The per- 
centage of low-mass stars that are fou nd to be active peak s at spectral type M5, where up to 
60% of stars may show Ha emission (IHawley et al.lll996l ). We also obtained high-resolution 
imaging of the field around SDSS-MEB-1 on 15 October, 2006 using the SCAM camera on 
Keck/NIRSPEC. In 0.36" seeing we find no companions between 1" and 25" down to an 
approximate limiting magnitude of i^=18.5. 



3. Analysis 

An initial inspection of the PAIRITEL photometry revealed large changes in the ob- 
served brightness of this system, a clear signature of an eclipsing binary . We estimated the 



period of the system using the Phase Dispersion Minimization method of lStellingweril (119781 ) 
in two stages. A first pass using coarse resolution in trial frequency indicated a prominent 
peak in the periodogram at around 0.4d with lesser peaks at double and half this period. 
We re-calculated a periodogram around P = 0.4d using a fine spacing in trial frequency of 
2.5 day -1 . We fit a parabola to the peak in this high-resolution periodogram to estimate the 
period of the system and found P=0.407037±0.000014d. We determined a precise time of 
zero phase with a parabolic fit to the center of the primary eclipse in the folded light curve 
and estimated the epoch of primary eclipse to be 2453988.7993 ± 0.0006 (HJD). 

We determined the radial velocities of both components of the binary by modeling the 
two Ha emission lines seen in the R ~ 3500 spectra. For both the R ~ 3500 and R ~ 900 
spectra we began by fitting Gaussian profiles to isolated emission lines in comparison lamp 
spectra in order to determine the width of the instrumental profile. We found that the 
individual Ha emission lines in the spectra of SDSS-MEB-1 were resolved in the R ~ 3500 
spectra so the Gaussian profiles were broadened to fit the data. We assume the same width 
for each of the Ha emission lines. Two of our R ~ 3500 spectra were acquired when 
the velocity separation between the two stars was close to maximal. In these spectra, the 
two emission lines are individually resolved and well separated and we fit a model of two 
super-imposed Gaussians to a 30A portion of the spectra around Ha. This model has five 
parameters: one for the centers of each Gaussian, one for the width of both Gaussians, one 
for their relative heights, and one for overall flux scaling. We used \ 2 minimization to fit for 
these parameters in the spectra taken near phase 0.25. A spectrum of SDSS-MEB-1 around 
the region of the two emission lines, along with the best-fit model, is shown in Figure [2j 
We used the fits to these well-separated lines to fix the relative heights of the emission lines 
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and their widths. The two emission lines are partially or totally blended in the rest of our 
data, particularly in the R ~ 900 spectra where they are both unresolved and blended. For 
the remaining R ~ 3500 spectra we used the fixed widths and relative heights determined 
previously to fit for the velocities of the two components. For the R ~ 900 spectra we used 
the relative heights determined previously and widths corresponding to the instrumental 
profile determined from fits to emission lines in comparison lamp spectra. We derived errors 
on the radial velocity estimates of each component in each spectrum based on the formal 
errors of the x 2 fit added in quadrature to the error estimate for the wavelength solution. 
All of the radial velocities are corrected to the solar system barycenter using a the IRAF 
routine bcvcorr. Our radial-velocity measurements are presented in Table 3 along with the 
estimated radial velocity curve parameters in Table 4. 

While the uncertainty on our individual radial velocity measurements, particularly those 
from the R ~ 900 spectra, are large, the radial velocities align very well with the epoch and 
period determined from the photometry alone. In order to estimate the system orbital 
parameters and the masses of each component we fixed the period and epoch of primary 
eclipse, assumed an eccentricity of e = 0.0, and fit a three parameter model to the radial 
velocity data. We estimated the two radial velocity amplitudes, K\ and K 2 , and the systemic 
velocity, V 1 with a \ 2 minimization. Assuming a Keplerian orbit with no eccentricity, the 
parameters of the radial ve locity orbit translate easily into estimates for M\ sin 3 i, M 2 sin 3 i, 



and a sin i ( iHilditchl 120011 ) . We estimated the errors on each of these values through a 
bootstrap simulation. We follow the convention that the primary star is the one eclipsed at 
phase 0, which is also the more massive star. This results in the marginally deeper eclipse 
depth seen at phase 0. The radial-velocity measurements, along with the best-fit model, 
are shown in Figure [31 Close inspection of the spectra resulting in two discrepant (off by 
~ 4<T from the model) measurements from UT Nov 21 (near phase 0.45) indicated that the 
measurements were likely biased by the superposition of a sky emission line on the shoulder 
of the marginally-blended Ha lines. In addition, some of our radial-velocity observations 
occur very close the time of primary eclipse. Assuming a rotation period locked to the 
orbital period of the system, we estimate that the amplitude of the Rossiter-McLaughlin 
effect is small comparable to our radial-velocity precision. 

In order to increase the power of light curve modeling, and to break a partial degen- 
eracy that exists in the determination of the ratio of the two stellar radii for systems with 
comparable eclipse depths or partial eclipses, we estimated the light ratio from our R ~ 3500 
spectra. We fit spectra taken near phase 0.0, when the component velociti es are perpendic 



ular t o our line of sight, to template M dwarf spectra made available by iBochanski et al. 



( 120071 ) . We found good agreement between our spectra and the active M4 dwarf template 



spectra. In Figure @] we show the spectrum of SDSS-MEB-1 compared to the template of 
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Bochanski et al.l (120071 ). We solved for the light ratio by fitting spectra taken at other phases 



to the sum of two active M4 dwarf template spectra Doppler shifted relative to each other 
and each scaled by a different multiplicative factor. The Doppler shifts of the two template 
spectra were fixed to values defined by the fit to the radial velocity curve. We calculated a 
X 2 goodness-of-fit between the model comprised of the weighted sum of the templates and 
our spectra for a range of relative brightnesses between 0.5 and 2.0. During the fitting, the 
area around the Ha emission lines was given no weight. A fit to the region around the 
minimum of the res ulting \ 2 surface yielded a light ratio estimate of L 2 /L\ = 0.85 ± 0.07. 



Baraffe et al.l (119981 ) presents theoretical evolutionary models of low-mass stars, the predic- 



tions of which we can compare to this estimate of the light ratio. The spectra used in our 
analysis span the wavel ength range 6 4 00A t o 7540A, comparable to the R band. At an age of 
2 Gyr and [M/H]=0.0, iBaraffe et all Jl998h predicts a light ratio of 0.744 (0.740 at t=5Gyr) 
in the Cousins/Bessel R band for the stellar masses determined from the fit to the radial 
velocity curve. Based on the reasonable agreement between our light r atio estimate and the 
value predicted from theory, we adopt the NIR light ratio estimates of IBaraffe et al.l (119981 ) 
throughout our modeling of the light curves. We estimated the error on this light-ratio esti- 
mate using a Monte-Carlo simulation by varying the masses (assuming sini = 1) of the two 
components of the binary within their errors. We estimated an uncertainty of ~ 0.2 in our 
light ratio estimat e and do not include a systematic contribution due to any biases that may 
be inherent to the IBaraffe et al.l (119 98) models themselves. 



We modeled our NIR light curves using the EB OP dPopper fc Etzell 1 1 98 lh technique as 
implemented in the code jktebop by J. Southworth (jSouthworth et al.l 120041 ) . Here, we fixed 
the mass ratio, period, eccentricity, epoch, and the brightness ratio L2/L1 and fit for the 
sum of the radii of the stars in terms of the orbital separation, (R\ + R2)/a 1 the ratio of the 
radii of the stars, the inclination angle i, and the ratio of the central surfaces brightnesses, 
J2I J\. We fixed the gravity darkening exponent to 0.32 as appropriate for a convective 
star (lLucylll967l ) and calculated the reflection coefficients from bolometric th eory. We used 
logarithmic limb darkening coefficients for a T e ff=3200K, log g = 5.0 star from IClaretl (120001 ) 
for NIR passbands. The jktebop implementation of EBOP allows for the estimation of errors 
on the parameters using Monte-Carlo simulations. Based on the quality and quantity of 
the photometric data collected with PAIRITEL, as well as the long time baseline of the 
observations, we did not attempt to model spots on the components of the SDSS-MEB-1 
system. Based on the EBOP models of the system, the out-of-eclipse variations are expected 
to be very small compared to the photometric errors. This allows us to evaluate the validity 
of our estimated photometric errors. We scaled our photometric error estimates until the 
reduced \ 2 °f the data at phase 0.7 < <fi < 0.9, assuming constant brightness, was 1.0. 
We derived the errors on our estimated parameters from 10 4 bootstrapping simulations and 
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include our estimate of the uncertainty on the orbital separation a. We found other sources 
of error from fixed parameters, such as the mass ratio and the limb darkening model, to be 
negligible. We fit the J—, H-, and K s — band light curves with consistent results. The 
photometric errors on our K s photometry are large (<Jk s > 0.15) compared to the eclipse 
depth and so the constraints on the system parameters from this band are relatively poor. 
We derived final values for the light curve parameters in two stages. First, the geometric (i, 
Ri + i? 2 , R2/R1) parameters were determined from weighted average of those from the J, 
H, and K s observations. With these geometric parameters fixed, we ran the EBOP analysis 
a second time to better determine the radiative parameters (central surface brightness ratio 
Jij J\ and total light ratio L2/L1). The model fit to our J-band data is shown in Figures 
[5] and [6] and the parameters from our EBOP analysis are shown in Table 5. The physical 
parameters of the SDSS-MEB-1 system derived from the radial velocity and photometric 
analyses are given in Table 6. We estimate the ratio of the radius of each component to the 
radius of its Roche lobe to be 0.389 and 0.367 for the primary and secondary, respectively. 

We estimated the effective temperatures (Teg) of each component of SDSS-MEB-1 by 
combining the results of the EBOP analysis with the Barnes- Evans Color-Surface Brightness 
relation. The EBOP analysis of the NIR light curves p rovides J2I J1 , the r atio of the central 
surface brightnesses of the two stars. From the work of iBeuermannl (120061 ) we can relate the 
ratio of the average surface brightnesses of the two stars to the ratio of their V — K colors. 
Since the temperatures of the two stars are expected to be similar, our assumption of the 
same limb darkening model for each star allows us to assume that the ratio of the central 
surface brightnesses is the sa me as the ratio of the a verage surface brightnesses. Following the 
color-temperature relations of iLeggett et al.l (119961 ) we convert NIR colors into temperatures. 
These relations are used to both turn the color ratio determined from the Barnes-Evans 
relation into a temperature ratio and to estimate the average temperature of the system. 
From the combined color of the SDSS-MEB-1, we estimate the luminosity-we ighted average 



tempe rature of the system. Converting the DENIS in to a Cousins I r followin g I Cost ado et al. 
(120051 ) and the 2MASS K s and J to K cit and J c u following ICarpenterl (120011 ) we estimate the 
luminosity- weighted average temperature of the SDSS-MEB-1 system to be 3330 ± 116f^. 
The temperature ratio estimate and the estimate of the average temperature of the system 
comprise a system two equations in two unknowns, 7\ and T 2 , which we solve numerically 
to arrive at the T e g estimates given in Table 6. The errors on the temperature estimates are 
derived from Monte Carlo simulations in which the errors from each of the colors and color 
transformations are propagated. We find that the temperatures of the two stars are consistent 
within the large errors and that the values agree well with the estimat ed temperatures for 
the c omponents of the system CM Draconis, which have similar masses ( jChabrier fc Baraffe 
1995h . 
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Discussion 



We combine the results of the light-curve and radial-velocity analyses to produce the 
physical parameters of the system listed in Table 6. The components of SDSS-MEB-1 fall 
in a mass regime were there are currently few direct mass and radius measurements from 
D EBs. Our measur e ment s are in good agreement with the empirical Mass-Radius relation 
of iBavless fc Oroszl d2006h as well as t he theoretical models (log* =8.0, 9.0; [M/H]=0) of 
Baraffe et al.l (119981 ). Recent work by IChabrier et al.l (120071 ) has explored the importance 
of rotation and magnetic activity in the evolution of low-mass stars in binaries. Since both 
components of SDSS-MEB-1 show Ha emission, demonstrating magnetic activity, a more 
precise set of system parameters could place important constraints on these new models. 
Our two new mass and radius measurements, along with models and other low-mass stars 
from the literature, are shown in Figure [71 If our system were very y oung then ou r radiu s 
estimates would would fall about 20% below the 0.1 Gyr relation of iBaraffe et al.l (119981 ). 
While SDSS-MEB-1 does lie in the vicinity of Orion, it is located more than 10° from 
known star-formin g regions. We also co mpare the log g and T e fj estimates from SDSS-MEB- 
1 to models from IBaraffe et al.l (119981 ). In Figure [8] we compare the two components of 
SDSS-MEB-1 to isochrones and mass-tracks. Using the solar metallicity models, we find a 
best-fit age for the syste m of logt > 8. 1 . We find good agreement between the mass-tracks 
and isochrones from the IBaraffe et al.l (119981 ) models. Both bias in the color-temperature 
relation we use and interstellar r eddening could be im portant in temperature estimates of 
DEBs. According to the maps of ISchlegel et al.l (119981 ). the total reddening in the direction 
of SDSS-MEB-1 in the NIR color we use to calibrate temperature is E(I — J) = 0.075. 
Based on comparison of the observed J-K colors of SDSS-MEB-1 to the same isochrones, 
we estimate a J— band modulus and a distance of the system of ~ 375 pc. At the position 
of SDSS-MEB-1 relative to the Galactic disk, we expect a significant portion of the total 
estimated reddening along the line of sight. Since the average color of the system is an 
important factor in the determination of the individual temperatures, for E(I — J) > we 
will underestimate the temperatures of the components of SDSS-MEB-1. We estimate that 
for each 0.1 mag increase in E(I — J) our estimated temperatures of the stars decreases by 
85K. 



With current data we lack sufficient precision to identify shortcomings in the models 
of low-mass stars. Increasing the total number of known low-mass DEB systems across a 
wide range of masses will improve our understanding of these objects. Deep, large-area 
surveys with sparse time sampling represent an excellent resource for the identification of 
new low-mass eclipsing systems. We have demonstrated that even as few as 15 high-quality, 
multi-color observations can be useful for identifying important new objects. While the 
amplitude of the eclipse signal is large enough to reliably identify promising candidates with 
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synoptic surveys, objects discovered in this way present unique difficulties in terms of follow- 
up observations. The sparse sampling of the discovery light curves requires a significant 
investment of observing time to generate a complete light curve. Focusing on relatively faint 
targets with good photometric precision results in a vast increase in the total number of 
targets, particularly M dwarfs, and affords an increase in the number of viable targets for 
searching for eclipsing binaries and transiting planets. At the same time, the spectroscopic 
and photometric followup becomes more problematic and more expensive as the candidates 
become fainter. The study of eclipsing systems with future synoptic surveys will require a 
shift in emphasis from discovery to follow-up observations. Robotic observing systems, such 
as PAIRITEL, are ideally suited for this type or work and will become an integral part of 
transient science produced by future synoptic surveys. 



5. Future Prospects 



The discovery of SDSS-MEB-1 highlights the potential importance of the SDSS-II, and 
future synoptic data sets, for studies of eclipsing systems. The SDSS-II data are qualita- 
tively different from the data sets typically used to search for low-mass eclipsing systems. 
Current planetary transit surveys tend to produce high-quality, high-cadence photometry for 
a relatively small number of targets. Future synoptic surveys will acquire fewer observations 
per target, with comparatively worse precision, but will observe orders of magnitude more 
targets. 

We attempt to quantify the potential impact of SDSS-II on studies of low-mass stars 
by estimating the number of DEB systems that are likely to be prese nt in the SDSS-I I 
data. We used the i - z to i-band absolute magnitude relations from IWest et al.l (120051 ) 
and color relations of iDavenport et al.l (120061 ) to estimate the absolute magnitudes of all of 
the point sources with i — z > 0.37 contained in the SDSS-II repeat scan region. These 
absolute magnitudes were first converted to Bessel ma gnitudes and then to mass estimates 
following the solar-age, solar- metallicity isochrones of iBaraffe et al.l (119981 ). Our estimate 
inclu des stars with mass 0.15Mfn < M < O.45M . We assume a total binary fraction of 
30% (IDelfosse et al.l I2004T) and integrate the mass-ratio distribution and semi- major axis 
distribution of iDuquennoy fc Mayorl (119911 ) in order to estimate the total number of binaries 
in our survey and the distribution of their parameters. We place several constraints on these 
binaries to estimate the number of systems that could be identified in the SDSS-II data 
and confirmed with follow-up observations. Since radial velocity observations are crucial 
for obtaining masses and radii of these stars, we constrain our targets to be brighter than 
i = 19.0 with both components exhibiting radial-velocity semi-amplitudes of at least 30 km 
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s -1 . In order to ensure that the system is double-lined, we also require that the ratio of the 
luminosities of the components is no less than 0.1. We require that the system be in eclipse 
for more than 10% of the time, so that the expected number of points in eclipse is ~ 2 with 
SDSS-II, and that the amplitude of the eclipse is ^ > 10%. 



We estimate the total number of stars in systems meeting these criteria in mass bins of 
width 0.05 M after scaling for the geometric probability of observing an eclipse for a given 
system. The results of our simulation are shown in Figure [9J We estimate that the SDSS- 
II may contain as many as 12 low-mass eclipsing binary systems where both components 
are in the mass range O.15M < M < O.45M . Our estimate attempts to count systems 
that are likely to be identified in the SDSS-II data and for which precise determination of 
mass and radius will be possible. Our simulation predicts the detection of between 1 and 7 
stars (90% confidence) of mass O.25M < M < O.3M , not inconsistent with our detection 
of the binary SDSS-MEB-1 (comprised of two low-mass stars). These estimates are only 
rough approximations since the binary f raction of M dwarfs, and the properties of their 
companions, are not well known. Work by lDelfosse et al.l (120041 ) points to an overabundance 
of short-period, q ~ 1.0 M dwarf binary systems. Any such systematic differences in the 
properties of M dwarf and G dwarf binaries would effect the validity of our estimates. 

The high quality, and shear quantity, of the data produced by synoptic surveys like 
SDSS-II also opens up the possibility of searching for small companions extending down 
into the regime of brown dwarfs and extrasolar planets. In principle, it is possible to detect 
eclipses o f Jupiter- mass o bject s in the SDSS-II data. The giant planet and brown dwarf 
models of iBurrows et al.l (120011 ) indicate that a companion more massive than lMj results 



in an eclipse depth of ^ > 10% for primaries with radii < 0.3R Q . While the 10% eclipse 
signal is probably at the limit of de tectability with our current photometric calibration of 
the SDSS-II data, llvezic et al.l (120071 ) have shown that homogeneous, 2% photometry can be 
derived from repeat SDSS observations. Future work focusing on purely differential photom- 
etry of stars in the SDSS-II data could potentially extend sensitivity to Neptune-size planets 
orbiting the smalle st stars. Extrapolating from th e 5 Gyr, so l ar-me tallic ity isochrones of 



Baraffe et al.l (119981 ) and the low- mass star colors of I West et al.l (120051 ) and iDavenport et al. 



j2006h (converting SDSS to Bessell colors) we estimate that our current selection criteria 
identifies ~ 40,000 {i < 19.0 ) targets with R± < O.3i? . While detection of a transit by 
a Jupiter-mass object may be feasible, follow-up radial-velocity observations to confirm the 
mass of the system will be very difficult for all but the brightest objects. Our experience 
with SDSS-MEB-1 indicates that radial velocities with 10 km s _1 precision should be feasi- 
ble with reasonable exposure times for targets brighter than i ~ 19.0. If we assume that a 
radial-velocity signal of semi-amplitude K = 30 km s" 1 can be detected with a reasonable 
investment of telescope time, then companions with M sin i > 95Mj can be detected in orbits 
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with periods P < 3 days around all of these targets. iBlake et al.l (120071 ) have demonstrated 
radial-velocity precision of 0.3 km s _1 using high-resolution infrared spectroscopy, but only 
for relatively bright objects (K < 12.0). Unfortunately, only about 1100 of the M dwarfs we 
selected are this bright but for these targets we are capable of radial- velocity confirmation of 
companions as small as Msini ~ 3Mj. Still, radial velocity followup of these fainter systems 
may be feasible with the next generation of 30m ground-based optical telescopes. 

Future deep, large-area synoptic surveys will contain a treasure trove of eclipsing systems 
and allow for the direct measurement of masses and radii of many low-mass objects down 
to the bottom of the main sequence. We demonstrate that photometry of faint objects with 
sparse time sampling can be used to discover important new eclipsing systems. With these 
systems, the challenge will lie in the follow-up observations necessary to produce precise 
estimates of the physical parameters of the system and it is likely that robotic observing 
facilities will play a crucial role in these investigations. We have presented here a new low- 
mass double lined eclipsing binary discovered in a data set that can serve as a prototype 
for the next generation of synoptic surveys, and demonstrated that follow-up observations 
of these objects are feasible. By extending to fainter objects and focusing on smaller stars 
it will be possible to greatly expand the number of viable targets for a search for eclipsing 
low-mass star and brown dwarf companions. 
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Fig. 1. — Multicolor light curve of SDSS-MEB-1 from the publicly available data in the 
first three seasons (2004 - 2006) of the SDSS II survey. Shown from top to bottom are the 
light curves in z-, i-, r-, and g-bands. The 2005 observations during two eclipses are readily 
apparent in all bands. The zero-points of the images were calibrated by tying the catalog 
PSF_COUNTS fluxes to SDSS DR4 catalog magnitudes. 
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Fig. 2. — Portion of the spectrum of SDSS-MEB-1 around the Ha emission lines after 
continuum subtraction. The best-fit model is overplotted as a dashed line. This R ~ 3500 
spectrum was taken at HJD 2454060.978 
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Fig. 3. — Keck radial- velocity measurements of the two components of SDSS-MEB-1. The 
radial velocities are derived following the methods outlined in Section 3. For this fit, the pe- 
riod and epoch of primary eclipse estimated from the PAIRITEL photometry were assumed. 
A barycentric correction has been applied to these radial velocities. The squares correspond 
to the radial velocities of Star 2 and the circles correspond to the radial velocities of Star 1 
(Mi > M 2 ). Star 2 eclipses Star 1 at phase 0.0 resulting in the slightly deeper eclipse depth 
seen at that phase in Figure 6. 




Fig. 4.— Keck LRIS (R~ 900) Spect rum of SDSS-MEB-1 compared to the M4 template 
spectrum from iBochanski et al.l (120071 ). Some broad molecular absorption features are indi- 
cated, along with atomic emission and absorption features (down arrows). 
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Fig. 5.— Top: PAIRITEL J-band differential light curve of SDSS-MEB-1. The best-fit 
model from an EBOP analysis of the light curve is overplotted. Bottom: Residuals of the 
model fit to the data. 
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Fig. 6. — Close up of the PAIRITEL J-band differential light curve around the two eclipses. 
The best-fit model from an EBOP analysis of the light curve is overplotted. The dashed 
line at J=15.87 is to guide the eye. The two eclipses are very similar in depth. Our best fit 
model is for a primary eclipse that is 0.015 mag deeper than the secondary. 
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Fig. 7. — Known low-mass stars with direct mass and radius measurements. The two 
co mponents of SDSS-MEB - 1 are p lotted as black s quare s. Other obje c ts tak en from Table 10 
oflLopez-Morales fc Riba3j2005h . lVaccaro et all J2007h . Eebb et all tood ). and the text of 



Lopez-Morales et al. (120061 ) are shown as gray circles. Overplotted are the emp irical relation 



of 



Bavless Sz Orosz (l2006h (short-dashed line), a 5 Gyr, [M/H]=0 mode l from iBaraffe et all 



(Il998f ) (solid line), a 0.1 Gyr [M/Hl=0 mo del from lBaraffe et all (Il998f ) (dotted line) and a 
5 Gyr, [M/H]=0 model from birardi et all ()2Q0J) (dash-dot line). 
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Fig. 8. — Th e comp onents of SDSS-MEB-1 compared to mass-tracks and isochrones from 



Baraffe et al.l (119981 ) models. The dashed lines are mass-tracks of masses 0.3, 0.25, and 0.2 



Mq. The solid lines are [M/H]=0 isochrones of ages log £=8.0, 8.3, and 9.0 years. From the 
combined light curve and radial velocity analyses, we estimate the masses of these two stars 
to be Mi = 0.272 ± 0.020M© and M 2 = 0.240 ± 0.022M© 
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Fig. 9. — Estimate of the total number of 0.15M Q < M < O.45M stars in DEB systems 
that can be identified in the SDSS-II repeat scan region and for which precise mass and 
radius determinations will be possible. The details of this simulation and the constraints 
placed on the targets are described in Section 4. Error bars correspond to 90% confidence 
interval assuming Poisson statistics. 
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Table 1. Basic Observables 



Parameter 


Value 


RA (J2000) 


03:18:23.88 


Dec (J2000) 


-01:00:18.4 


9 


20.818±0.041 


r 


19.290±0.015 


i 


17.681±0.007 


z 


16.792±0.010 


J 


15.396±0.048 


H 


15.725±0.056 


K s 


14.623±0.101 


id 


16.907±0.100 



Note. — Position and bright- 
ness of SDSS-MEB-1. The 
J,H, and K s values are from 
2MASS and the SDSS val- 
ues are from t he SDSS cat- 
alog llAdelman-McCarthv et al.l 
2006). The i]j magnitude is 
from the DENIS survey. 



Table 2. NIR Photometry 



HJD-2400000 


Band 


Mag 


17 


53995.8909 


.1 


15.350 


0.043 


53995.8929 


J 


15.356 


0.045 


53995.8950 


J 


15.350 


0.043 


53995.8971 


J 


15.288 


0.043 


53995.8992 


J 


15.422 


0.067 


53995.9013 


.1 


15.295 


0.045 


53995.9034 


J 


15.323 


0.064 


53995.9056 


J 


15.393 


0.048 


53995.9077 


J 


15.406 


0.048 


53995.9098 


J 


15.543 


0.051 



Note. — PAIRITEL NIR photometry of 
SDSS-MEB-1. Errors reported here have 
been scaled as described in Section 3. Only 
a portion of the photometry is reported here. 
All of the data are available in the on-line 
edition. 
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Table 3. Radial-velocity Measurements 



HJD-2400000 


RVi 


°\RVl 


RV 2 


a RV 2 


setup 




km s _1 


km s~ x 


km s _1 


km s -1 




54000.13646 


15.79 


21.70 


-159.32 


10.33 


400/8500 


54000.14036 


27.63 


24.28 


-120.50 


11.14 


400/8500 


54000.14427 


14.60 


29.89 


-137.36 


13.29 


400/8500 


54060.84514 


-49.19 


14.90 


-64.03 


7.39 


1200/7500 


54060.84904 


-55.76 


19.44 


-65.21 


8.83 


1200/7500 


54060.97794 


-147.97 


5.88 


49.48 


5.19 


1200/7500 


54060.98185 


-152.38 


5.79 


45.47 


5.18 


1200/7500 


54061.02872 


-117.77 


7.52 


-20.83 


5.50 


1200/7500 


54061.03263 


-114.74 


7.67 


-28.96 


5.52 


1200/7500 



Note. — Radial-velocity measurements of the two components of SDSS- 
MEB-1. These values were estimated following the procedures outlined in 
Section 3. The last column indicates the groves mm (400 or 1200) and the 
blaze angle (7500A or 8500A) of the grating used. 



Table 4. Radial Velocity Fit Parameters 



Parameter 



Value 



V 7 (km s- 1 ) 

Ki (km s- 1 ) 

K 2 (km s- 1 ) 

O— C rms (km s _1 ) 

e (fixed) 



-61.24±2.31 
107.51±5.22 
121.70±4.22 
13.96 
0.0 



Note. - System parameters for 
SDSS-MEB-1 estimated from modeling 
spectroscopic observations as described 
in Sections 2 and 3. 
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Table 5. Light-Curve Fit Results 



Parameter 


J Band 


H Band 


K s Band 


Adopted Value 


(ill + «2)/a 

Ra/Ri 

Ri/a 

R 2 /a 

i(°) 


0.279±0.006 
0.928±0.019 
0.144±0.0034 
0.134±0.0034 
85.87±0.59 


0.284±0.008 
0.911±0.022 
0.149±0.0050 
0.135±0.0039 
85.56±0.35 


0.270±0.022 
0.870±0.045 
0.140±0.012 
0.126±0.012 
85.72dz0.73 


0.280±0.005 
0.916±0.014 
0.145±0.003 
0.134±0.003 
85.65±0.28 


Ja/Ji 

L2/L1 (input) 
ia/^l (best-fit) 
LD [c,d] (fixed) 


0.977±0.029 
0.792±0.230 
0.842 
[0.5619,0.5212] 


0.980±0.048 
0.811±0.200 
0.814 
[0.5471,0.5310] 


1.075±0.100 
0.818±0.200 
0.813 
[0.4595,0.4603] 




P (day) 
T min (HJD) 
O-C rms (mag) 


0.06 


0.08 


0.14 


0.407037±0.000014 
2453988.7993±0.0006 



Note. — System parameters for SDSS-MEB-1 estimated from modelin g of the photo metric data. 
Limb darkening (LD [c,d]) parameters are for the logarithmic model of IClaretl l|2000l ). The value 
of J2 is the relative central surface brightness of the secondary (Ji = 1.0). T m i n is the time when 
the primary is eclipsed. The O — C values indicate the rms of the residuals of the data minus the 
model. 
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Table 6. Physical Parameters for SDSS-MEB-1 



Parameter Value 



Mi (M Q ) 


0.272±0.020 


M 2 (M Q ) 


0.240±0.022 


a(R e ) 


1.850±0.047 


Rt (R ) 


0.268±0.0090 


R 2 (R©) 


0.248±0.0084 


log 9i 


5.016±0.043 


log 92 


5.029±0.050 


T X 


3320 ± 130 K 


T> 


3300 ± 130 K 


M V i 


11.55 ± 0.30 


My 2 


11.86 ± 0.37 



Note. — System param- 
eters for SDSS-MEB-1 esti- 
mated from modeling of the 
photometric and spectroscopic 
followup described in Sections 
2 and 3. The values for sur- 
face gravity, logg, are given in 
cgs units. The estimates of the 
absolute \/-band magnitudes of 
the stars are based on fits to 
iBaraffe et all lll998h models 



